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Abstract 

Let A{t) denote the cluster produced by internal diffusion limited 
aggregation (internal DLA) with t particles in dimension d > 3. We 
show that A{t) is approximately spherical, up to an 0{\/\ogt) error. 

In the process known as internal diffusion limited aggregation (internal 
DLA) one constructs for each integer time t > an occupied set A{t) C Z'^ 
as follows: begin with ^(0) = and ^(1) = {0}. Then, for each integer 
t > 1, form A(t + 1) by adding to A(t) the first point at which a simple 
random walk from the origin hits Z'^ \ A{t). Let Br C R'^ denote the ball of 
radius r centered at 0, and write B,. := H Z'^. Let io^ be the volume of 
the unit ball in M"^. Our main result is the following. 

Theorem 1. Fix an integer d > 3. For each 7 there exists an a = a{'y,d) < 
00 such that for all sufficiently large r, 

P {b,_,^ C Ai^dT^) C B,+,^}' < r-\ 

We treated the case d = 2 in [.JLSIO] (see also the overview in [.JLS09]), 
where we obtained a similar statement with logr in place of \/logr. To- 
gether with a Borel-Cantelli argument, this in particular implies the follow- 
ing [JLSIO]: 

Corollary 2. The maximal distance from dB^ to a point in one (hut not 
both) of o-nd A{ijJdr'^) is a.s. O(logr) when d = 2 and O(yTogr) when 
d>2. 
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These results show that internal DLA in dimensions d > 3 is extremely 
close to a perfect sphere: when the cluster A(t) has the same size as a ball 
of radius r, its fluctuations around that ball are confined to the \/log r scale 
(versus logr in dimension 2). 

In [JLSIO] we explained that our method for d = 2 would also apply in 
dimensions d >3 with the log r replaced by \/log r. We outlined the changes 
needed in higher dimensions (stating that the full proof would follow in this 
paper) and included a key step: Lemma A, which bounds the probability of 
"thin tentacles" in the internal DLA cluster in all dimensions. The purpose 
of this note is to carry out the adaptation of the d = 2 argument of [JLSIO] 
to higher dimensions. We remark that in [.JLSIU] we used an estimate from 
[LBG92] to start this iteration, while here we have modified the argument 
slightly so that this a priori estimate is no longer required. 

One way for A{u!di^'^) to deviate from the radius r sphere is for it to 
have a single "tentacle" extending beyond the sphere. The thin tentacle 
estimate [.JLSIO, Lemma A] essentially says that in dimensions d > 3, the 
probability that there is a tentacle of length m and volume less than a small 
constant times m'^ (near a given location) is at most e~'^"^ . By summing 
over all locations, one may use this to show that the length of the longest 
"thin tentacle" produced before time t is 0{\^logt). To complete the proof 
of Theorem 1, we will have to show that other types of deviations from the 
radius r sphere are also unlikely. 

Lemma A of [-JLSIO] was also proved for d = 2, albeit with e~^™ replaced 

g-cm /iogm_ However, when d = 2 there appear to be other more "global" 
fluctuations that swamp those produced by individual tentacles. (Lideed, 
we expect, but did not prove, that the logr fluctuation bound is tight when 
d = 2.) We bound these other fluctuations in higher dimensions via the 
same scheme introduced in [.JLS09, -JLSIO], which involves constructing and 
estimating certain martingales related to the growth of A(t). It turns out 
the quadratic variations of these martingales are, with high probability, of 
order logt when d = 2 and of constant order when d > 3, closely paralleling 
what one obtains for the discrete Gaussian free field (as outlined in more 
detail in [JLSIO]). The connection to the Gaussian free field is made more 
explicit in [JLSll]. 

Section 1 proves Theorem 1 by iteratively applying higher dimensional 
analogues of the two main lemmas of [JLSIO]. The lemmas themselves are 
proved in Section 3, which is the heart of the argument. Section 2 contains 
preliminary estimates about random walks that are used in Section 3. 
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A brief history of internal DLA fluctuation bounds 

The history of fluctuation bounds such as the one in Corollary 2 is as follows. 
In 1991, Lawler, Bramson, and Griffeath proved that the limit shape of inter- 
nal DLA from a point is the ball in all dimensions [LBG92]. In 1995 Lawler 
gave a more quantitative proof, showing that the fluctuations of A^u^r'^) 
from the ball of radius r are at most of order 0(r^/^ log'' r) [Law95]. In 
December 2009, the present authors announced the bound O(logr) on fluc- 
tuations in dimension d = 2 [JLS09] and gave an overview of the argument, 
making clear that the details remained to be written. In April 2010, As- 
selah and Gaudilliere [AG 10a] gave a proof, using different methods from 
[JLS09], of the bound ©(r^/^^^+i)) in all dimensions, improving the Lawler 
bound for all d > 3. In September 2010, Asselah and Gaudilliere improved 
this to 0((logr)^) in all dimensions d>2 with an O(logr) bound on "inner" 
errors [AG 10b]. In October 2010 the present authors proved the O(logr) 
bounds (announced in December 2009) for dimension d = 2 and outlined the 
proof of the 0{y/logr) bound for dimensions d > 3 [JLSlOj. In November 
2010, Asselah and Gaudilliere gave a second proof of the 0{\/logr) bound 
[AGlOc]. Their proof uses methods from [AGlOl)] along with Lemma A of 
[JLSIO] to bound "outer" errors and a new large deviation bound (in some 
sense symmetric to Lemma A) to bound "inner" errors. 

More references and a more general discussion of internal DLA history 
appear in [JLSIO]. 

1 Proof of Theorem 1 

Let m and i be positive real numbers. We say that a; G Z'^ is m-early if 

X £ A{ujdi\x\ - m)'^), 

where ojd is the volume of the unit ball in M*^. Likewise, we say that x is 
Mate if 

x^AM\x\+£y). 
Let <fm[7'] be the event that some point of A{T) is m-early. Let Ci[T] be 
the event that some point of ^(^r/uiay/'^-i ^® ^-late. These events correspond 
to "outer" and "inner" deviations of A(T) from circularity. 

Lemma 3. (Early points imply late points) Fix a dimension d > 3. For 
each 7 > 1, there is a constant Cq = Co{j,d), such that for all sufficiently 
large T, if m > Co\/TogT and i < m/Co, then 

p(£:„[r]nA[T]^) <r-io^. 
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Lemma 4. (Late points imply early points) Fix a dimension d > 3. For 
each 7 > 1, there is a constant Ci = Ci{'y,d) such that for all sufficiently 
large T, if m > £ > Cl^/^ogT and i > Ci{{logT)my/^ , then 




I eh 4 = c(riogr)i/3 



Figure 1: Let be the smallest m' for which A{T) contains an m' early 
point. Let l^ be the largest i' for which some point of B^j,^^^y/d_^/ is /-late. 
By Lemma 3, {i'^,m'^) is unlikely to belong to the semi-infinite rectangle 
in the left figure if ^ < m/Co. By Lemma 4, {i'^,m'^) is unlikely to belong 
to the semi-infinite rectangle in the second figure if ^ > C7i((logr)m)V3. 
Theorem 1 will follow because > mo = T is impossible and the other 
rectangles on the right are all (by Lemmas 3 and 4) unlikely. 



We now proceed to derive Theorem 1 from Lemmas 3 and 4. The lemmas 
themselves will be proved in Section 3. Let C = max(Co, Ci). We start with 

mo = T. 

Note that A{T) C By, so F{£t[T]) = 0. Next, for j > we let 
£j = max(C7((logr)mj)V3^C70ogr) 

and 

TTlj+i = C£j. 

By induction on j, we find 

p(f™^.[r]) <2jT-io^ 
P(£.,[T]) < (2j + l)r-io^ 
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To estimate the size of £j ,letK = log T and note that < , where 
Then 

< max(i^l/3+V9+■■■+l/3^^l/3^i^V2) 
so choosing J = log T we have 

^1/3^ < 2 

and 

Ij < 2K^/^ < C^J\ogT. 

The probability that A{T) has £j-late points or mj-early points is at 
most 

(4J + l)T-i°^ < T-^^' < r~^. 

Setting T = ujd^'^, (■ = (-j and ra = mj, we conclude that if a is suffi- 
ciently large, then 

which completes the proof of Theorem 1 . 

2 Green function estimates on the grid 

This section assembles several Green function estimates that we need to 
prove Lemmas 3 and 4. The reader who prefers to proceed to the heart of 
the argument may skip this section on a first read and refer to the lemma 
statements as necessary. Fix d >3 and consider the d-dimensional grid 

G = {(xi, . . . , Xd) G M'^ : at most one Xj ^ Z}. 

In many of the estimates below, we will assume that a positive integer k and 
a y G Z"^ have been fixed. We write s = \y\ and 

n = n{y,k) ■.= gnBs+k\{y}- 

For X G 90, let 

P(X) = Py,k{x) 

be the probability that a Brownian motion on the grid Q (defined in the 
obvious way; see [JLSIO]) starting at x reaches y before exiting Bs+k- Note 
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that P is grid harmonic in (i.e., P is linear on each segment of \ Z*^, 
and for each x G 17 n Z'^, the sum of the slopes of P on the 2d directed edge 
segments starting at x is zero). Boundary conditions are given by P{y) = 1 
and P{x) = for X E [dO) \ {y}. The point y plays the role that C, played 
in [JLSIO], and P plays the role of the discrete harmonic function H(^. One 
difference from [ ' ] is that we will take y inside the ball (i.e., A; > 1) 
instead of on the boundary. 

To estimate P we use the discrete Green function g{x), defined as the 
expected number of visits to x by a simple random walk started at the origin 
in Z'^. The well-known asymptotic estimate for g is [Uch98] 



g{x) - ad\x\'^ 



<C\x\-'^ (1) 



for dimensional constants and C (i.e., constants depending only on the 
dimension d). We extend 5 to a function, also denoted g, defined on the grid 
Q by making g linear on each segment between lattice points. Note that g 
is grid harmonic on ^ \ {0}. 

Throughout we use C to denote a large positive dimensional constant, 
and c to denote a small positive dimensional constant, whose values may 
change from line to line. 

Lemma 5. There is a dimensional constant C such that 

(a) P{x) < C/(l + |x-y|^-2)_ 

(b) P{x) < Ck{s + k + l- \x\)/\x - y\'^, for \x-y\> k/2. 

(c) max P{x) < Ck/{s - r - k)'^^^ for r < s - 2k. 

Proof. The maximum principle (for grid harmonic functions) implies Cg{x — 
y) > P{x) on ri, which gives part (a). 

The maximum principle also implies that for x G (7, 

P{x)<C{g{x-y)-g{x-y*)) (2) 

where y* is the one of the lattice points nearest to {s + 2k + Ci)y/s. Indeed, 
both sides are grid harmonic on fi, and the right side is positive on dBgj^k 
by (1), so it suffices to take C = {g{0) — g{y — y*))~^ . 
Combining (1) and (2) yields the bound 

Ck 

P(x) < TTTT) foi' l^; — y| > 2k. 

\x — y\"-~^ 
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Next, let z G dBg+k be such that 1^ — y| = 2L, with L > 2k. The bound 
above imphes 

Ck 

P{x) < for X G Bl{z) 

Let z* be one of the lattice points nearest to {s + k + L + Ci)z/\z\. Then 

F{x) = adL^-'' -g{x-z*) 

is comparable to L^"*^ on dB2L{z*) and positive outside the ball Bl{z*) (for 
a large enough dimensional constant Ci — in fact, we can also do this with 
Ci = 1 with L large enough). It follows that 

P{x) < C{k/L^-^){L'^-'^)F{x) 

on d{B2L{z*) n ri) and hence by the maximum principle on B2l{z*) H il. 
Moreover, 

F{x) <C{s + k + l- \x\)/L'^-^ 
for x a multiple of z and s + k — L < \x\ < s + k. Thus for these values of x, 

P{x) < C{k/L)F{x) < Ck{s + k + l- \x\)/L'^ 

We have just confirmed the bound of part (b) for points x collinear with 
and z, but z was essentially arbitrary. To cover the cases \x — y\ < 2k one 
has to use exterior tangent balls of radius, say k/2, but actually the upper 
bound in part (a) will suffice for us in the range |x — y| < Ck. 

Part (c) of the lemma follows from part (b). □ 

The mean value property (as typically stated for continuum harmonic 
functions) holds only approximately for discrete harmonic functions. There 
are two choices for where to put the approximation: one can show that 
the average of a discrete harmonic function u over the discrete ball is 
approximately u(0), or one can find an approximation Wr to the discrete 
ball Bf. such that averaging u with respect to Wr yields exactly u{0). The 
divisible sandpile model of [LP09] accomplishes the latter. In particular, the 
following discrete mean value property follows from Theorem 1.3 of [LP09]. 

Lemma 6. (Exact mean value property on an approximate ball) For each 
real number r > 0, there is a function Wr '■ l^"^ ^ [0, 1] such that 

• Wr{x) = 1 for all X £ ^r-c, for a constant c depending only on d. 

• Wr{x) = for all x ^ Bj.. 
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For any function u that is discrete harmonic on B^, 

Wr{x){u{x) -n(0)) = 0. 



The next lemma bounds sums of P over discrete spherical shells and 
discrete balls. Recall that s = \y\. 

Lemma 7. There is a dimensional constant C such that 
(a) ^ P{x) < Ck for allr < s + k. 

xeBr+i\Br 



(b) 



(c) 



j;(P(x)-P(0)) 



< Ck for all r < s. 



{p{x)-pm 



xeB 



s + k 



< Ck\ 



Proof. Part (a) follows from Lemma 5: Take the worst shell, when r = s. 
Then the lattice points with |x — y| < k, s < \x\ < s + 1 are bounded by 
Lemma 5(a) 

s'-'^s'^-'ds = k 

Jo 

(volume element on disk with thickness 1 and radius k in Z'^"^ is s'^~'^ds.) For 
the remaining portion of the shell, Lemma 5(b) has numerator k{s + k — s) = 
k"^, so that 

poo 

/ k'^s-'^s'^-^ds = k 

Jk 

Next, for part (b), let Wr be as in Lemma 6. Since P is discrete harmonic 
in Bs, we have for r < s 

J2 Wr{x){P{x) - P{0)) =0. 

Since Wr equals the indicator 1b,. except on the annulus B^ \ B,._c) and 
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\wr\ < 1) we obtain 

5](P(:r)-P(0)) 



< Yl \Mx)\\p{x)-pm 

< Ck. 

In the last step we have used part (a) to bound the first term; the second 
term is bounded by Lemma 5(b), which says that P{0) < Ck/s'^~^. 

Part (c) foUows by spUtting the sum over Bg+fc into k sums over spherical 
shells Bg+j \ '^s+j-i for j = 1, . . . , /c, each bounded by part (a), plus a sum 
over the ball B^, bounded by part (b). □ 

Fix a > 0, and consider the level set 

U = {x eG\ g{x) > a}. 

For X G dU, let p{x) be the probability that a Brownian motion started at 
the origin in G first exits U at x. 

Lemma 8. Choose a so that dU does not intersect Z*^. For each x G dU , 
the quantity p{x) equals the directional derivative of g/2d along the directed 
edge in U starting at x. 

Proof. We use a discrete form of the divergence theorem 

divF = Vi^f; • (3) 
^ du 

where y is a vector-valued function on the grid, and the integral on the left 
is a one-dimensional integral over the grid. The dot product vjj-V is defined 
as Cj ■ V{x — Ocj), where Cj is the unit vector pointing toward x along the 
unique incident edge in U . To define the divergence, ioi z = x + tej, where 
<t <l and X eZ'^, let 

d 

diyV{z) := —Cj ■ V{z) + (5^(z) ^(e^- • V{x + Oej) - Cj ■ V{x - Oej)). 

J i=i 

If / is a continuous function on U that is on each connected compo- 
nent olU — TJ^, then the gradient of / is the vector- valued function 

v = vf = {df/dxi,df/dx2, . . . , df/dxd) 
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with the convention that the entry df / dxj is if the segment is not pointing 
in the direction Xj. Note that V/ may be discontinuous at points of Z*^. 

Let G = —g/2d, so that div VG = Sq. If u is grid harmonic on U, then 
div Vtt = and 

div {uVG - GVu) = u{0)6o. 

Indeed, on each segment this is the same as {uG' — u'G)' = u'G' — u'G' + 
uG" — u"G = because u and G are hnear on segments. At lattice points u 
and G are continuous, so the divergence operation commutes with the factors 
u and G and gives exactly one nonzero delta term, the one indicated. 

Let u{y) be the probability that Brownian motion on U started at y first 
exits U at x. Since u is grid-harmonic on U, we have divVu = on 
hence by the divergence theorem 

u{Q)= I diw{uVG-GVu)=^uvu-yG. □ 



Next we establish some lower bounds for P. 
Lemma 9. There is a dimensional constant c > such that 

(a) P(0) > ck/s'^-^. 

(h) Let k = I, and z = {1 - ^)y. Then 

min P{x) > c/m'^~^. 

a;£B(z,m) 

Proof. By the maximum principle, there is a dimensional constant c > 
such that 

P{x) > c{g{x -y)- ad{k/2f-'') 
for X G -6^/2(2/)- III particular, 

P{x) > ck^-'^ for all \x - y\ < k/A 

Now consider the region 

U = {x^g: g{x) > aais'f-"} 

where s' is chosen so that \s' — {s — k/%)\ < 1/2 and all of the boundary 
points of U are non-lattice points. (A generic value of s' in the given range 
will suffice.) 
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By (1), this set is within unit distance of the bah of radius s — k/8. Let 
p{z) represent the probabihty that a Brownian motion on the grid starting 
from the origin first exits U at z £ dU. Thus 

n(0) = u{z)p{z) (4) 

for all grid harmonic functions u in U. 

Take any boundary point of z G dU. Take the nearest lattice point z* . 
Let Zj be a coordinate of z largest in absolute value. Then \ zj\ > \z\/d. The 
rate of change of in the jth direction near z has size > l/(i|2;|'^~"^, 

which is much larger than the error term Cl^l"*^ in (1). It follows that 
on the segment in that direction, where the function g{x) — ad{s — k/8)'^~'^ 
changes sign, its derivative is bounded below by l/2d\z\'^~^ . In other words, 
by Lemma 8, within distance 2 of every boundary point of z G dU there 
is a point z' E dU for which p{z') > c/s'^~^. There are at least ck'^~^ such 
points in the ball Bfc/4(y) where the lower bound for P was ck'^~'^, so 

P(0) > ck'^-'^k'^-^/s'^-^ = ck/s'^-\ 

Next, the argument for Lemma 9(b) is nearly the same. We are only 
interested in A; = 1. It is obvious that for points x within constant distance 
of y (and unit distance from the boundary at radius s + 1, the values of P{x) 
are bounded below by a positive constant. We then bound P{{s — 2m)y/\y\) 
from below using the same argument as above, but with Green's function 
for a ball of radius comparable to m. Finally, Harnack's inequality says 
that the values of P{x) for x in the whole ball of size m around this point 
(s — 2m)y/\y\ are comparable. □ 

3 Proofs of main lemmas 

The proofs in this section make use of the martingale 

M{t)=My,,it):= iPix)-Pm 
xeAy,fc(t) 

where Ay^i.{t) is the modified internal DLA cluster in which particles are 
stopped if they exit 0. As in [JLSIO], we view Ay^k{t) as a multiset: points 
on the boundary of where many stopped particles accumulate are counted 
with multiplicity in the sum defining M. In addition to these stopped par- 
ticles, the set Ay^k(t) contains one more point, the location of the currently 
active particle performing Brownian motion on the grid Q. 
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Recall that P = Py k and M = My^k depend on A;, which is the distance 
from y to the boundary of $7. We will choose k = \ for the proof of Lemma 3, 
and k = ai for a small constant a in the proof of Lemma 4. Taking /c > 1 is 
one of the main differences from the argument in [JLSIO]. 

Proof of Lemma 3. The proof follows the same method as [JLSIO, Lemma 
12]. We highlight here the changes needed in dimensions d> 3. We use the 
discrete harmonic function P{x) with k = 1. Fix z G Z'^, let r = \z\ and 
y = {r + 2m)z/r. Let 

Ti = \ujd{r - 

where o;^ is the volume of the unit ball in W^. \i z is m-early, then z € A{Ti); 
in particuar, this means that r > m, so that r+m, r+2m are all comparable 
to r. Since A; = 1, we have by Lemmas 5(c) and 9(a) 

P(0) ^ l/r^-\ 

where « denotes equivalence up to a constant factor depending only on d. 
First we control the quadratic variation 

N 

S{t) = lim Y.^M{t.{) - M(t,_i))2 

on the event that there are no (m + l)-early points by time T. 

As in [JLSIO, Lemma 9], there are independent standard Brownian motions 
i?^, . . . such that each increment [S(n + 1) — S{n))\^^^^\^Y is bounded 
above by the first exit time of B"^ from the interval [— a^, 6n]) where 

an = P(0) « ^ 

hn = max P(x) < — ; — —-r, , . 

\x\<(n/uaV/d+m+i [r + 2m- {{n/LOdfl'^ + m+ l)^-^ 

Here we have used Lemma 5(b) in the bound on 6„. 

Unlike in dimension 2, we will use the large deviation bound for Brownian 
exit times [JLSIO, Lemma 5] with A = cm? instead of A = 1. Here c is a 
constant depending only on d. Note that bn < l/m°'~^, for all n < Ti, so 
this is a valid choice of A in all dimensions d > 3 (that is, the hypothesis 
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logE 



^(a„ + bn) < 3 of [JLSIO, Lemma 5] holds). We obtain 
r 1 

n=l 

(r + m - {n/oJdY/'^ - "^Y^^ ^ 
r'^ ^ {r + m — J — 1)'^ ^ 



< 



I [r + m — j — 1)" ^ 



Note that the last step uses d > 3. Taking A = cm for small enough c we 
obtain 



E 



Therefore, by Markov's inequality, 

p({5(ri) > i/c} n £m+i[Tf ) < e"^-'"' < (5) 

Fix z G By and t G {1,...,T}, and let Qz^t be the event that z G 
A{t) \ A(t — 1) and z is m-early and no point of A(t — 1) is m-early. This 
event is empty unless (t/Wf^)^/'^ + m < < (t/w^)^/'^ + m + 1; in particular, 
the first inequality implies t <Ti. We will bound from below the martingale 
M{t) on the event Qz,t H Ce[TY. With no £-late point, the ball ^r-m~e-i 
is entirely filled by time t. Lemma 7(b) shows that the sites in this ball 
contribute at most a constant to M{t) (recall that k = 1). The thin tentacle 
estimate [JLSIO, Lemma A] says that except for an event of probability 
g-cm ^ iiiQ^Q aj-g order m'^ sites in A{t) within the ball B(z, m). By Lemma 
9(b), P is bounded below by c/m'^~^ on this ball, so these sites taken together 
contribute order m to M{t). Each of the remaining terms in the sum defining 
M{t) is bounded below by — P(0), and there are at most ir'^~^ sites in 
A{t) \ ^r-m-i-i- So these terms contribute at least 

-ir'^-^{l/r'^-^) = -£> -m/C 

which cannot overcome the order m term. Thus 



n {M^{t) < m/C} n CeitT) < e-'"" • (6) 
We conclude that 

p(Q,,t n Ci[TY) < P(Q,,t n {S{t) > i/c}) 

n {M{t) < m/C} n CeW) 
+ mS{t) < l/c}n{M{t) > m/C}). 
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The first two terms are bounded by (5) and (6). Since M{t) = B{S{t)) for 
a standard Brownian motion B, the final term is bounded by 

P \ sup B{s) > m/C \ < e-'=W^)V2 < ^-20^^ □ 

lo<s<l/c J 



Proof of Lemma 4- Fix y G Z'^, and let L[y] be the event that y is £-late. 
Let s = \y\, and set A; = a£ in the definition of P. Here a > is a small 
dimensional constant chosen below. Note that the hypotheses on m and 
i imply that i is at least of order \/log T; after choosing a, we take the 
constant Ci appearing in the statement of the lemma large enough so that 
> 10007logr. 

Case 1. 1 < s < 2fc. Then P{0) ^ l/s'^~^. Let 

With a„ = -P(O) and bn = 1, we have S(n + 1) — 'S'('^) < ^n, where r„ is the 
first exit time of the Brownian motion B^ from the interval [— a„, (Note 
that because we take = I, the indicator 1£-^^i[t]= is not needed here as it 
was in the proof of Lemma 3.) We obtain 

Ti 

logEe'^('^i) < ^logEe^" < riP(O). 

t=i 

Let Q = riP(O). By Markov's inequality, P(5(ri) > 2Q) < e'^ . 

On the event L[y], the site y is still not occupied at time Ti. Accordingly, 
the largest M(Ti) can be is if Ay^k(Ti) fills the whole ball B^+fc (except for 
y), and then the rest of the particles will have to collect on the boundary 
where P is zero. The contribution from ^s+k is at most Ck"^ by Lemma 7(c). 
The number of particles stopped on the boundary is at least 

Ti-2u;dis + kf> y. 

Therefore, on the event L[y] we have 

M{Ti) < Ck^ - Y^'(O). (7) 

Note that Q := TiP(O) « {s + 1)'^ / s'^'^ > i'^/{k/2Y-^, so by taking a = kjl 
sufficiently small, we can ensure that the right side of (7) is at most —Q/A. 
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Also, Q>f> 10007logT. Since M{t) = B{S{t)) for a standard Brownian 
motion B, we conclude that 



F{L[y]) < F{S{T,) > 2Q) + ¥ ^^mf^^B{s) < -Q/4 
< e-'3 + e-(0/4)V4Q 



Case 2. s > 2k. Then by Lemma 5(c) with r = 1, and Lemma 9(a), we 
have P(0) ^ k/s'^'^. First take 

To = [ujd{s + k- 2,mY\ 

(or To = if s + A; — 3m < 0). As in the previous lemma (but taking A = 1 
instead of A = cm?) we have 



logE 



/o {s + k-{n/ojdYl'^) 
The last inequality follows from d > 3 and m > k/a. By Markov's inequality, 

p({5(ro) > c + A;2} n f,^[r]=) < e-^'' < 

Now since 

(Ti - ro)P(O) ^ ms'^-\k/s'^-^) = km 

we have 

logEe^(^i)"^(^o) < Cfem. 

Thus (since km > k"^) 

¥[{S{Ti) > 2Ckm} n £m[TT) < 2T-'^°^ . (8) 

As in case 1, the martingale M{Ti) is largest if the ball B^+fc is com- 
pletely filled, and in that case the total contribution of sites in this ball is 
at most Ck'^. On the event L[y], the number of particles stopped on the 
boundary of at time Ti is at least 



Ti - #B,+,, > ojdiis + i)''-is + k + C)'^) « is 



d-l 
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Each such particle contributes — P(0) ~ —k/s'^ ^ to M(Ti), for a total 
contribution of order —k£ = —k'^/a. Taking a sufficiently small we obtain 
M{Ti) < Ck'^ - k'^/a < -k'^. We conclude that 

F{L[y] n Smim < n{S{Ti) > 2Ckm] n £^[m+ 

+ ¥{{S(Ti) < 2Ckm} n {M(Ti) < -k'^}). 

The first term is bounded above by (8), and the second term is bounded 
above by 

p| inf i?(s)<-/fc2\ <e-fcV4Cfcm^^-207_ 
ls<2Ckm J 

Hence F{L[y] n£m[TY) < ST-^Ot. Since Ce[T] is the union of the events L[y] 
foi y & B := B^j.^^^ji/<i_£, summing over y £ B completes the proof. □ 
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